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■ First Hints of Jet Quenching at RHIC 

O 

^ ; A. Drees, SUNY at Stony Brook 

At this conference first data from RHIC has been presented. Spectra of charged hadrons 
and identified neutral pions obtained in central collisions exhibit a depletion at large trans- 
verse momenta compared to expectations deduced from pp and pp data and lower energy 
. heavy ion data. While spectra measured in peripheral collisions exhibit the expected 
power-law shape, spectra from central collisions are closer to exponential. In addition, a 
significant azimuthal anisotropy of high momentum charged particle production has been 
Q\ ; found. All observations are in qualitative agreement with theoretical predictions that 
^ ■ quark matter formed in heavy ion collisions quenches jet production. 

1. Introduction 

o : 

Parton-parton scattering with large momentum transfer, so-called hard scattering pro- 
cesses, provide unique probes for studing heavy ion collisions. The reason is simple: hard 
i— h \ scattering occurs early in the collision, well before quark matter is expected to form. The 
scattered partons will sense the full space-time evolution of the collision volume and thus 
Ch ■ probe the later formed hot and dense phase. 

>• . Scattering with large momentum transfer either results in the production of high mass 
particles, like charm quarks of which a small faction bind in J/ip charmonium states, or 
in high momentum quarks or gluons which fragment into jets of hadrons. If the scattered 
partons penetrate quark matter significant modifications of J/i/j production and jet frag- 
mentation are expected. In fact, the discovery of charmonium suppression U is one of 
the corner stones of the argument that quark matter has been formed already at CERN 
energies @. Jet production should also be altered significantly. High momentum partons 
should lose a significant fraction of their energy by gluon bremsstrahlung effectively sup- 
pressing jet production. This phenomenon, commonly refered to as "jet quenching", was 
predicted a decade ago ||. 

Among the large number of produced particles present in the final state of a heavy ion 
collision jets can not be reconstructed directly. However, one of the jet fragments will 
always carry a major fraction of the jet momentum. These so-called leading particles 
manifest themselves in a power-law shape of the transverese momentum distribution. If 
jets are quenched a depletion of the high momentum tail of the spectra is expected. In 
addition, due to the binary nature of hard scattering processes, most jets are produced 
in pairs, thus azimuthal correlations between high p± particles might also serve as an 
experimental observable. 

At CERN SPS energies measurements of high p± particle production where performed 
H||. A conclusive interpretation failed so far, mostly because the data can be explained 
by hard scattering || as well as transverse flow [0. In contrast, first data from RHIC show 
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characteristic features consistent with the anticipated jet quenching. This talk summarizes 
the experimental evidence shown at QM2001. In the next two sections data from nulceon- 
nucleon collisions are evaluated and used to establish a reference to which RHIC data can 
be compared. Section 4 compares inclusive charged particle production to the reference. 
Additional information on identified particle production, and its implicitions are discussed 
in section 5. The final section summarizes the talk. 



2. Elementary reactions 

Searching for new phenomena in heavy ion collisions requires a detailed understanding 
of elementary nucleon-nucleon collisions. Data on high momentum particle production 
from pp or pp collisions exists for various beam energies, but not for the initial RHIC 
energy of \fs nn = 130 GeV. Fig.[l] shows inclusive charged particle p± spectra obtained 
from pp collisions at the ISR ||, and from pp collisions by UAl || at CERN, and by CDF 



at FNAL |10[|. With increasing beam energy high pj_ particle production is enhanced, 



reflecting the increase of the jet production cross-section. 
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Figure 1 . Inclusive charged particle production 
for pp and pp collisions [^|4jo| . The ISR data for 
7T, K, p, and p were added to obtain the charged 
spectrum. All data were fitted with a power- 
law function (dotted lines) and interpolated to 
t/s = 130 GeV (thick line), systematic errors of 
the interpolation are shown as dashed lines. 
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Figure 2. Charged particle and neutral pion 
data from different experiments [^11 1^] at y/s 
from 500 to 630 GeV. The consistency of the data 
is used to estimate the systematic errors on the 
absolute normalization of individual experiments. 



To generate a reference p± distribution for the initial RHIC energy y/s nn = 130 GeV, the 
data shown in Fig.|l| were fitted with the empirical functional form d 2 a jdp\ = Aj (po+p±) n - 
The cross sections were then interpolated to 130 GeV at several fixed p±. Finally, the 
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interpolated cross sections were fitted with the identical power-law function to obtain a 
smooth reference distribution. The fit parameters obtained are A = 2.75 10 5 , po = 1.71, 
and n = 12.42. The fits and the interpolation are shown in the figure. The dashed 
lines below and above the interpolation indicate the systematic uncertainty which results 
from (i) uncertainties in the interpolation procedure which increases with pj_ and more 
importantly (ii) systematic discrepancies of data sets take by different experiments. The 
latter error was estimated to by about 20% by comparing different data sets at similar 
y/s. An example is given Fig. ^| |pHTl|]. 

The same magnitude of systematic errors on the nucleon-nucleon reference is estimated 
by comparing to other extrapolations done by the STAR collaboration |~2],[l3| and X.N. 



Wang WM as well as comparing it to theoretical calculations of the cross section [|I5|,|T6 



3. Extrapolating to AA and the Nuclear Modification Factor 

To extrapolate from pp to AA collisions non-trivial assumptions need to be made. Lets 
first assume all high p±_ particle production results from binary hard collisions. In the 
absence of nuclear effects the cross section should then simply scale by the atomic number 
squared A 2 . For a specific impact parameter selection corresponding to a fraction / of 
the inelastic nucleus-nucleus cross section the average number of binary collisions 

of the event sample (Nonary) can be calculated from the nuclear overlap integral and the 
inelastic nucleon-nucleon cross section aT el 



nn 
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In the following values of a™^ 1 = 42 mb and a AA l = 7200 mb are used consistently. 
Experimentally, the number of binary collisions can be evaluated from measured global 
observables like the forward energy or the event multiplicity by a Glauber model approach 
(e.g. in [18|]). To search for deviations from the nn-reference one may divide thep_i_ spectra 
from heavy ion collisions by the scaled nn-reference. This ratio, denoted as Raa, will be 
referred to as the nuclear modification factor following the suggestion from [III]. 



D / v da AA dyd 2 p ± 

Raa{P±) = ttt rj -n , 2 (2) 

(A 'binary) d(T pp / dyd Z p± 

If simple scaling with the number of binary collisions holds true the nuclear modification 
factor R A a should be unity. At low momenta, say below 1 GeV/c, Raa must be smaller 
than one since in this p± region the cross section is expected to scale with the number of 
participating nucleon rather than with the number of binary collisions. Thus scaling with 
{Nunary) is expected only at high p±, say above 2-3 GeV/c. 

Below p± of 10 GeV/c the situation is complicated by known nuclear effect. Already 
in the late 70's it was discovered by Cronin et al. that high p± particle production in 
p-nucleus collisions is enhanced beyond simple binary scaling. Traditionally this enhance- 
ment, now called Cronin effect, has been parameterized as a p A = A a( - P± ^ a pp . Fig.^ gives a 



compilation of the p± dependent exponent from different fixed target experiments [19,20 



Though, so far there is no quantitative theoretical description of the Cronin effect, it is 
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Figure 3. Compilation of data for the Cronin expo- 
nent a(p±) |l!||2(J. The excess above binary scaling, 
the "anomalous nuclear" enhancement, sets in around 
2 GeV/c when a(p±) increases above unity. 



Figure 4. Nuclear modification factor for 
charged and pion data obtained in 158 AGeV Pb- 
Pb and Pb-Au collisions at CERN HHIH (figure 
5 from [0). As shown in 14 the nuclear modifi- 
cation factor is consistent with the Cronin effect. 



commonly accepted that it originates from initial state multiple scattering, i.e. multiple 
scattering of parton in the collision ETl . 

Data taken at 



The Cronin effect has also been observed in heavy ion collisions [23 



CERN has been recently compiled and analyzed in terms of the nuclear modification 




factor [|T|]. FigJ4| reproduces a figure from |14| showing Raa for Pb-Pb collisions. Raa 
continuously increases, crosses unity around 1.5 GeV/c and eventually saturates above 
2.5 GeV at a value Raa ~ 2. At higher energies a similar behavior is expected, though 



one might expect a reduction with increasing energy 

In summary, if previously observed phenomena are included, at RHIC Raa should 
be above the naive binary scaling, which probably underestimates Raa and below the 
empirically observed value from CERN, which might overestimate Raa- This conclusion 
implicitly assumes that (i) the relative particle abundances - in particular at high p± - 
do not change going from nucleon-nucleon to nucleus-nucleus extracted from RHIC data 
collisions, and that (ii) collective radial flow as well as (iii) shadowing do not alter the 
spectra at high p± significantly. 



4. Inclusive charged particle data from RHIC 

Fig.|5] compares the inclusive p± distribution of negatively charged hadrons from STAR 
[O] to the spectra of all charged hadrons presented by PHENIX |2l 



Both data sets are 

obtained from the 5% most central Au-Au collisions at pseudo-rapidity zero. The data are 
consistent over the entire range observed out to 5 GeV/c. This is remarkable, given the 
preliminary nature of the data and the known difficulties of absolute measurements. There 
are subtle difference between the data sets which amount less than 40% and are largest 
around 2 GeV/c. This difference is used to estimate the systematic error on the absolute 
normalization. In the following the error is assumed to be ±20%, which is consistent with 
the errors quoted by the experiments. 

To compare to the nn-reference the nuclear modification factor is calculated. The nn- 
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Figure 5. Comparison of data on negatively 
charged hadrons obtained by STAR Jl| and 
charged hadrons from PHENIX Q for the 
most central 5% of the collisions. Both data 
sets are independently absolutely normalized. 
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Figure 6. Nuclear modification factor determined 
as ratio of the data from figure ^| over the nn-reference 
described in section 2, scaled by (Nonary) = 945 (see 
text for details). The thin line indicates the estimated 
systematic uncertainty. Also given is the expectation 
for scaling with {Nbinary) (dashed line) and the nuclear 
modification factor deduced from CERN data (dash- 
dotted line) 



STAR used a 



reference was scaled by (N binary ) = 945 ± 140 determined by PHENIX [[ 
value of (Taa) — 26 mb -1 , corresponding to (N binary ) = 1092 |L3| to scale their reference. 
The difference of both values is used as systematic uncertainty of 15% on (N binary ). This 
estimate is consistent with the ~ 15% error given by PHENIX. 

The resulting nuclear modification factor Raa is shown in Fig||. On a linear scale the 
subtle differences between the data sets become more apparent. The systematic errors 
on Raa are indicated by the full line, they results from (i) the error on the nn-reference 
(~20% but p± dependent), (ii) a 20% error on the normalization, and (iii) a 15% error on 
(Nbinary)- The errors quoted here are somewhat larger than those evaluated by PHENIX 
and STAR, the reason is two fold: (i) the errors take into account the differences between 
the data and (ii) the error on the nn-reference was estimated more conservatively. 

Initially, Raa increases up to a of about 2 GeV/c where it saturates at a value 
of 0.6 to 0.8. At high the nuclear modification factor Raa seems to decrease again. 
In the region from 2 to 3 GeV Raa might be consistent with unity within the rather 
large systematic uncertainties. However, it is not converging towards one at higher p±, 
as suggested by simple binary scaling. The figure also indicates Raa observed at the 
CERN-SPS (dash-dotted line). Obviously, Raa at RHIC is significantly below the value 
found for CERN SPS data. 

PHENIX has also presented a systematic survey of different centrality selections. The 
data are shown in Fig.[7]. The spectra for the most peripheral collisions exhibit a pro- 
nounced power-law shape which seems to vanish for the most central selection. The 
significant difference between central and peripheral collisions is more clearly visible in 
the ratio, which is depicted in Fig.[8[ To compare to the binary scaling assumption and 
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Figure 7. Inclusive charged particle production measured by PHENIX j24|] for different centrality 
selections (specified in the figure). All data sets are absolutely normalized. The lowest data set represents 
the most peripheral selection and the highest data set the most central selection. 



to Raa the cross sections were normalized by the number of binary collisions for the spe- 
cific centrality selection. The large uncertainty in {Nunary) for the peripheral data sets 
(~ 60%) is reflected by the error band also shown in the figure; this ia an uncertainty on 
the absolute scale of the ratio but not on its p± dependence. 

Both ratios, central 5% to 80-92% and to 60-80% show the same trend: a rise up to p± of 
1.5-2 GeV/c followed by a decrease at higher momenta. Such a behavior is qualitatively 
in agreement with a power-law spectrum becoming more exponential with increasing 
centrality. Within the systematic errors the ratio could be close to one, but as in Fig.|| 
the ratio does not converge towards unity at large p± as suggested by the binary scaling 
scenario. If one assumes that peripheral collisions resemble nucleon-nucleon interactions 
one can interpret the ratio shown in Fig.|] as an alternative measure of the Raa- The 
ratio derived here and Raa shown in Fig.|6| agree reasonably well. Although the systematic 
errors have similar size they have very different origin. 

That something interesting is happening at high p± is corroborated by data on az- 
imuthal angular correlations presented by STAR [Z5]|. The data is presented in Fig 



in 



terms of the p± dependence of i^. Here V2 is the second Fourier coefficient of the azimuthal 
track density distribution measured with respect to the reaction plane. The coefficient 
t>2 measures an anisotropy of particle production relative to the reaction plane, which is 
typically interpreted as elliptic flow. At low p±_ vi continuously increases at a rate pre- 
dicted by model calculations based on hydrodynamics [p6[ . Around 2 GeV/c the pattern 
changes and t> 2 seems to saturate. At high p±, say above 2 GeV/c, particle production can 
no longer be described by hydrodynamics and consequently the anisotropy is expected 
to vanish. Interestingly, the saturation of v 2 and Raa occur at similar p± suggesting a 
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Figure 8. Ratio of central to peripheral data 
on charged particle production (as shown in Fig.0). 
The lines indicate the systematic error band, which 
results mostly from the uncertainty in the number 
of binary collisions for the peripheral event selec- 
tions (N binary ). 
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Figure 9. STAR data Q on azimuthal 

anisotropy — V2 — measured versus the reaction 
plane as a function of p±_ . Also shown are recent 
model calculations invoking a quark matter phase 
]27| . The three different lines indicate different ini- 
tial gluon densities. 



possible connection of both phenomena. 

Angular correlations between jets should not contribute to v 2 because jet production 
is uncorrelated to the reaction plane. If, however, jets loose significant energy in the 
dense medium an anisotropy will be introduced by the geometry of the reaction volume. 
In transverse direction, perpendicularly to the beam axis, the reaction volume will be 
almond shaped for non-zero impact parameter. Thus, the average path length of a jet 
in the dense medium will be smaller in plane than out of plane. It is plausible that jets 
produced out of plane will loose more energy than those produced in plane, effectively 
leading to an anisotropy with respect to the reaction plane at medium and high momenta. 

The figure compares the data to a recent theoretical calculation [[^j which combines a 
hydrodynamic calculation with perterbative QCD, and modeles jet production quenched 
by energy loss. The different curves correspond to different initial gluon densities; for a 
density of 500 the model simultaneously describes also the p± spectra reasonably well. 



5. Identified pion data from RHIC 

Additional information on high p±_ production of neutral pion, proton and anti-proton 
is available from PHENIX fl28"|--j30|j . These first data show that in Au-Au collisions a much 



smaller fraction of the high p± charged particles are pions than in pp collisions. Minimum 
bias data for 7r~, K~ and p is presented in Fig.[L0]. Above 2.0 GeV/c the ratio p/7T~ 
(as well as the p/ir + ratio) approaches unity, significantly larger than the value of ~ 0.2 
observed in pp collisions at the ISR [J In addition, the spectra suggest that pj ti~ continues 
to increase towards higher p±, quite different from m± scaling found in pp collisions. Large 
collective radial flow in heavy ion collision might offer a possible explanation [|3T ,32|. 



Independent of the interpretation, at high p± the particle compositions in Au-Au at 
RHIC deviates from pp collision where the ratio charged /it was about 1.6 above 1.5 GeV/c. 
From Fig.|irJ one deduces charged/ir of about 2 at 1.5 GeV/c, increasing to ~ 2.5 at about 
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Figure 10. Inclusive 7r~, K~ , and p pro- 
duction from minimum bias Au-Au collisions pre- 
sented by PHENIX [§§]. All three distributions 
are absolutely normalized. 
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Figure 11. Neutral pion transverse momentum 
distribution measured by PHENIX [§|]. The data 
is compared to the appropriately scaled nucleon- 
nucleon reference for pion production (systematic 
errors given as dashed lines). 



2.2 GeV/c. Thus Raa measured for charged particles should not be interpreted as nuclear 
modification factor for pions. Raa underestimates the depletion of pion production by 
approximately a factor ~1.5. 
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Figure 12. Ratio of central to peripheral 

neutral pion data. The thin lines indicate the 
systematic uncertainty of the data. 
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Figure 13. Comparison of neutral pion pro- 
duction to various model calculations from |l6[ ] 
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The measurement of neutral pions by PHENIX |28] reveals directly the more substantial 
suppression of high p± pion production. Two data sets are given in Fig.[TT], a 10% central 
and a peripheral (60-80%) event selection. The systematic uncertainty of the absolute 
normalization is included in the error bars. 

Both data sets are compared to the nucleon-nucleon reference. Here the reference, 
which was deduced for charged particles in section 2, was scaled down by the charged /ir 
ratio of 1.6 observed at the ISR. The systematic error bands shown in the figure contain 
the uncertainty of the reference (see section 2), a 10% error on the charged/ V scale 
factor, as well as 15% and 60% error on the number of binary collisions for central and 
peripheral data, respectively. While the peripheral data are well described by the reference 
distribution, the central data are significantly below. From the ratio of the data to the 
nn- reference one deduces an average nuclear modification factor of Raa ~ 0.4 for neutral 
pions in central collisions. This value is smaller by a factor of 1.5 compared to Raa 
obtained for charged particles, consistent with the large p and p yield at high p±. For the 
peripheral data, Raa is about unity reflecting the good agreement of the data and the 
wi-reference. 

As alternative measure of Raa gives the ratio of central to peripheral data. For 

all p± the ratio is below unity, indicating the substantial suppression of the yield in central 
collisions relative to the binary scaling assumption. Very similar to the charged particle 
data the ratio peaks between 1.5 and 2 GeV/c and decreases towards higher p±. 

Finally in figure 0, the neutral pion data are compared to theoretical model calculations 
by X.N. Wang JTH]. Three different calculations are shown: (i) an estimate of the scaled 
wi-reference, (ii) a scaled reference modified by modeling the Cronin effect (here called 
p^-broadening) and shadowing, as well as (iii) a calculation inferring "jet quenching" 
by introducing an average energy loss of 0.25 GeV/fm for the scattered partons. The 
peripheral data agrees reasonably well with all three scenarios while the central data are 
significantly below the calculations not including energy loss. 



6. Conclusion and Outlook 

Only a few months after the end of its first Au-Au run, RHIC has produced a surpris- 
ingly rich sample of data, reaching out in p± well into a region where a significant contri- 
bution of hard scattering processes is expected. In central collisions a depletion of high 
p±_ particle production is discovered. This depletion is observed in charged particle pro- 
duction measured by STAR and PHENIX as well as in neutral pion data from PHENIX. 
The same depletion is found when comparing central to peripheral data. Additional in- 
formation results from azimuthal anisotropy of particle production. All observations are 
consistent with a scenario in which quark matter formed during the collision suppresses 
jet production. 

Whether these first hints for "jet quenching" will hold true remains to be seen. Draw- 
ing definite and quantitative conclusions from the present data certainly is premature 
in particular in view of the large systematic uncertainties. Besides better control of the 
absolute normalization it will be essential to measure the nn-reference in the same exper- 
iments. Data at higher p± will help to disentangle soft physics like radial flow from hard 
scattering. 
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At the end I would like to express my gratitude to many who have helped preparing my 
summary talk - before and during this exciting conference. I specially thank T.Ullrich, 
P.Jacobs, and R.Snellings for providing the STAR data. X.N. Wang generously provided 
calculations and figures for my talk. Also many colleges in PHENIX deserve credit, in 
particular G.David, F.Messer, and J. Velkovska. 
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